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ABSTRACT 

This paper generalizes and extends the concept of the Schuler oscillation that occurs in the theory of 
inertial navigation systems, allowing one to see how the Schuler phenomenon affects inertial navigation 
systems operating in space. We show why a low earth orbit satellite's orbital period is identical to the 
period of the Schuler pendulum, which is the period of the errors for terrestrial inertial navigation systems. 
We also show that the generalized form of the Schuler oscillation takes the same fonn as the Hill-Clohessv- 
Wiltshire equations for satellite relative motion, and that the period of the out-of-plane motion in 
neighboring satellite relative trajectories is the same as the Schuler period. Finally, we describe how INS 
gyro drift manifests itself in different coordinate systems for the orbital case. These results may assist 
orbital flight dynamics and attitude control systems engineers in the design and analysis of INS-equipped 
spacecraft 


INTRODUCTION 

People working in the field of inertial navigation systems (INS) are familiar with the phenomenon of 
Schuler oscillations [1] in INS. The majority of them are familiar with the concept of the Schuler 
pendulum. However, the fact that the period of both is identical to the period of low earth orbit (LEO) 
satellites and to the period of the oscillatory- change of the relative location of two adjacent LEO satellites is 
not widely known in the INS community, if at all. It is also not widely known that when the INS is installed 
on a LEO satellite traversing a near-circular orbit, the INS error equations are the same as the Hill- 
Clohessy- Wiltshire equations that describe the change of the relative location of two adjacent LEO 
satellites. On the other hand, orbital mechanics professionals are well aw r are of the orbital period of LEO 
satellites and the equations of their relative motion, but in general they are unfamiliar with the Schuler 
pendulum and Schuler oscillations in INS, and thus are unaware of the identity between all four periods. In 
this paper we examine the four phenomena; namely, the Schuler pendulum, Schuler oscillations, the LEO 
satellite orbital period, and the period of the relative motion between two adjacent LEO satellites. In 
particular we show why the LEO satellite’s orbital period is identical to the well-known INS Schuler 
period. We also show- that if the INS error equations are generalized to a non-terrestrial case, a generalized 
form of the Schuler oscillation exists, which takes the form of the Hill-Clohessy-Wiltshire equations for 
satellite relative motion. 

SCHULER PENDULUM 

In a 1 923 paper [1], while claiming that INS could not be realized because “differentiation between the 
force of gravity and acceleration is impossible”, the German Physicist, Maximilian Schuler, explained that 
if one can build a pendulum that has the length of the earth’s radius then one could move at will the 
pendulum point of support and still have the pendulum point in the direction of the vertical. (Schuler 
explained that this stemmed from the fact that the center of mass of the apparatus remained at the center of 
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the earth and thus at rest.) It is very well known that. T p , the period of such mathematical pendulum, is 
given by 

T,-2^ (!) 

where R is the radius of the earth, and g is its gravity. 

SCHULER OSCILLATIONS IN INS ERROR PROPAGATION 

The simplest way to describe the generation of oscillations in INS is to examine the behavior of a 
stable-platform INS in response to tilt error. 

For simplicity, consider a stable-platform INS where the platform control system is programmed 
to keep the platform horizontal with one of the platform axes, say the x-axis, pointing north, the y-axis 
pointing east and the z-axis pointing in the direction of the local nadir. Consequently the preferred platform 
coordinate system is the North, East, and Down system. Suppose that there is a platfonn tilt error, denoted 
by cp , about the platform North-axis (see Fig. 1). As a result the East accelerometer reads gravity as 

acceleration component f = — g • sill tp (the sign is negative because the accelerometer reads gravity as 

acceleration in the opposite direction). The INS computer is set to integrate tins measured f to obtain 
velocity. Thus the computed East velocity is the integral of f; in other words. 



LOCAL 
EAST * 


v = -g*sincp 


(2 a) 


Usually (p - 0, therefore vve can write 


v = -gcp (2.b) 

In order to keep the platfonn level, the platform control system rotates the platform about the North axis by 
the angular rate V / R , thus 


V 



( 3 ) 


0 





Differentiation of the last equation and substitution of V from Eq. (2.b) into the resultant equation yields 

<P+|-<P = 0 (4) 

The last equation is that of a harmonic oscillator whose period, T ws , is given by 



(5) 


This period, known as Schuler period, is well known in Inertial Navigation technology. A comparison 
between Eqs. (5) and (1 ) reveals that 


Tins “ T p (6) 

It is only ironical that this period, a characteristic feature of INS, is named after Schuler who maintained 
that the construction of INS was “an entirely impossible undertaking” [1,2]. 

LEO SATELLITE ORBITAL PERIOD 

Consider a LEO satellite on a circular orbit at the surface of a hypothetical perfectly smooth 
airless earth. In order to be in orbit the satellite has to travel at the orbital velocity where the centrifugal 
acceleration balances gravity. That is, the orbital velocity, V , has to satisfy the relation 


II- 

R =S 

thus 

vfyiR 


(7) 

( 8 ) 


The distance covered after one orbit is the circumference of earth, which is 2~R . Therefore write: 


VT 0 = 2-R 


(9) 


T 0 being the duration of one revolution. Substituting V from Eq. (8) into the last expression yields 


which gives 


VgR-T 0 = 2*R 



( 10 ) 
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It means that the orbital period also equals the Schuler period! 

This can also be shown in the following way. It is well known that the Keplerian parameter, n , 
also known as Mean Motion, is given by 


n = 



( 12 ) 


but 


p- R u- 1 

” ” R 3 R 2 R R 2 R 


( 13 ) 


where R is the vector of earth radius from the center of earth to the satellite, and 1 R is a unit vector in the 
direction of this vector, therefore 


From Eqs. (12) and (14) we obtain 


and since n = 2% IT 0 , we obtain 


which is the Schuler period. 



(14) 

(15) 


(16) 


GENERALIZATION OF THE SCHULER OSCILLATION IN INS ERROR PROPAGATION 


The above derivation uses a simplified scalar INS error model to describe the Schuler oscillatioa 
This section generalizes that explanation to the vector case. It then shows how the generalized form of the 
INS error equations results in Hill’s, equations. 

The INS senses a specific force, f , consisting of the true acceleration, and the acceleration of 
gravity, g , where the latter is interpreted as a negative acceleration, thus: 


f=a-g(R t ) (17) 

where gravity is a vector function of the true position, R t . The INS accelerometer measurement, f m , is 
corrupted by bias (among other errors neglected here), 5 a , thus: 

f m =a-g(R t ) + 5a (18) 

As a result, the INS reconstruction of the true acceleration is as follows: 

a c =f ra + g(R 0 ) (19 a) 

where R c is the computed position vector. Using Eq. (1 8) the last equation becomes 


a 0 =a-g(R t ) + g(R c ) + Sa 

(19. b) 

a, = a + 5a + [g(R t + 5R)-g(R,)] 

(19 c) 


where 5R = R c - R t . Then using a truncated Taylor series expansion yields 
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a = a + 5a + 


(19. d) 


Consequently 


s(R,) + J! «R-s(R.) 


MR f 


a = a +5a + — | 
dR 


5R + ... 


Therefore, to compute its position, the INS performs the following integration: 

R c = jj a c dt dt = j J (a + 5a +G(R t )8R + ...) dt dt 


where G(R t ) = <5g/ 8 R| r . Since (see Bq. 13) 


pR 

R 5 ” 


p R 
R^R 


= -gl 


R 


where 1 R 


[0 


0 


1] T 


. it can be shown in a straightforward manner that 


In reality, the true position is 

R t = jj adtdt 

To get a differential equation for the INS errors, insert Eq. (23) into Eq. (20) to get 

R c = R t + jj(5a + G(R,)5R + ...)dtdt 

and differentiate Eq. (24) with respect to an inertial frame. This yields 

R 0 = R t + Sa +G(R t )5R + ... 

which, using the definition of 5R , can be written as: 

5R -G(R,)5R =5a + ... 

Substitution of G (R t j from Eq. (22) into the last equation yields: 


(19e) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


5R + |-(l-31 R l^)5R = 5a+... (27) 

Normally m terrestnal INS equations the 2 -a>:is of the reference coordinates points along the vertical, 
therefore 


:> 


and Eq. (27) becomes 


(l-31 R l*)*diag0 1 -2} 


(28) 


5R + — [1 1 -2]5R*5a (29) 

R 

The horizontal; that is, the x and y components of Eq. (29) form hannonic oscillators, with a natural 
frequency given by the. square root of the scalar coefficient of the gravity gradient, and a corresponding 

period = 2%*Jg~/ R , which is identical to the expression given in Eq. (5). 

THE PERIOD OF HILL'S EQUATION 

Consider a LEO satellite in the vicinity of another LEO satellite. We name the first satellite 
Primary Satellite ( PS) and the second Secondary Satellite (SIS}. Define an orbital coordinate system whose 
origin is in the PS. Its x-axis is in the direction of the satellite zenith, the v-axis is in the direction of the 
satellite velocity vector, and the 2 -axis is nonnal to the orbital plane on the side that completes the set to a 
right-hand coordinate system (see Fig. 2). The differential equations that describe the relative displacement 
of the SS with respect to the PS in this coordinate system are known as HilFs equations or Clohessy- 
Wiltshire equations [3]. In the absence of external forces these equations are: 


X 



Fig. 2: The Orbital Coordinate System for Relative LEO. Motion 


x-3o3 2 x-2oy = 0 

(30. a) 

y+2oi = 0 

(30. b) 

z + O) 2 z = 0 

(30.c) 


All three components have an oscillatory mode at the angular frequency CO (for the r- component this is 
obvious). It has been shown [3] that 
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(we approximate here the radius of the circular orbit by the earth radius) and as we saw before* T, the 
period that corresponds to such angular frequency is 


T = 2tt 



(32) 


In fact it is easy to see without mathematical developments that this is the period of the z component.. For 
simplicity assume that the orbit of the PS is circular and that the initial conditions on the x 9 v, and z 
components, and their derivatives are as follows X 0 = x o = y o = y o = 0 , z o = 0 * z o & 0 . Then the 

orbit of the SS is also a circular orbit that intersects the orbit of the PS twice in an orbit as shown in Fig. 3. 
Obviously the distance 


X 



Fig. 3: The Relative Location of the Secondary Satellite with Respect to 
the Primary Satellite along the z-axis. 


between the two satellites changes direction once on even' orbit and the frequency at which this happens is 
the orbital frequency which, in the preceding section, was found to be the Schuler frequency. 

EXTENSION OF THE INS ERROR EQUATIONS TO NON-TERRESTRIAL MECHANIZATION 

If we consider an INS installed on a satellite circling the earth in a near-circular orbit, and if 
Eq. (26) is written in a local-vertical local-horizontal frame, Hill’s Equations result. To see this, consider 
the following general INS position error equation [4]: 


rj2 


— T-5R + 2 
dt 2 


f r 1 > \ 

© r x— 5R 

dt 


+ -^-© r x 5R + ©' x (e> r x 5R) - G(R)5R = Aa - y x f (33 .a) 


where C0 r is the angular velocity of frame r vrith respect to an inertial frame, and the pre-superscripts of the 
d/dt operator indicate the frame in which the derivative is taken. The vectors Aa and V| / denote. 
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respectively, the accelerometer errors and the angular error caused by the gyro errors which are denoted by 
the vector 8 . (The latter errors are known as *psi 5 errors). The c psL angular error develop according to [4] 


— \j/ + <v> r XV | l = s 


(33 .b) 


At rest on earth C0 r contributes a 24 hour oscillatory mode in the propagation of q/ , which through the 
term -\j/ x f on the right side of Eq. (33. a) contributes this mode to the propagation of 5R . In our case of 
a LEO satellite G) r is the orbital angular velocity at which the satellite orbits the earth, and according to 
Eq. (II) the orbital period corresponding to 0 r is exactly the Schuler period. As a result, the oscillatory 
mode in the c psP angular error has a Schuler period rather than a 24 hour period. Because CO r is constant 
then the third term on the left of Eq. (33.a) drops, and so is the term -\|/ * f on the right side of Eq. (33. a) 
because in orbit f = 0 . Therefore, if we consider only accelerometer bias, Eq. (3 3. a) becomes: 


rj2 

—t-5R + 2 

v -- „ 

where, as before, 5 a , is a vector whose three components are the components of the three corresponding 
accelerometer biases. Let us denote the magnitude of (D r by n, that is ri = Jo> r | . In a circular orbit (see 

Eq. 1 5) n* = g / R , so Eq. (22) can be written as 

G(R)5R = -n 2 5R 4- 3ti 2 1 r 1 r SR (34) 


co T x — 5R 
dt 


+ G> r x ^o r x 5R^ - G(R)5R = 5a (33. c) 


Choose the axes of the rotating frame such that 0t) r = nk r , and the x- axis is pointing in the direction of 
the zenith; that is, 1 R = i r . Also, let 


5R = xi r + yj r + zfc r 


(35) 


Inserting Eqs. (33 .c), (34), and (35) into Eq. (26) results in 


x - 2ny - 3n 2 x = 5a x + ... 

(3 6. a) 

y + 2nx = 5a v + ... 

(36.b) 

z + n 2 z = 8a z +... 

(36.c) 


in which w r e recognize the homogenous part as the Hill-Clohess) r -Wiltshire equations of relative motion, 
forced by the accelerometer bias and higher-order terms. We conclude that when on a circular orbit, the 
INS position error due to accelerometer bias behaves like the relative position between two adjacent 
satellites on low earth orbit In particular^ rhe relative position along the normal to the PS orbital plane 
oscillates at the Schuler frequency. Moreover, the c psi 5 angular emor too has a Schuler period, which, 
unlike in terrestrial INS, is not coupled into the INS position error. 
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Remark: In some spacecraft the INS reference coordinate system is ah inertial one. In such cases equation 
(3 3. a) reduces to 1 d\|/ / dt = S where 1 designates the inertial frame. In such cases all three of the INS 
channels exhibit unbounded error growth due to gyro drift, and the spaceflight INS must use external aids 
such as star trackers to periodically re-align its attitude. 

CONCLUSIONS 

In this paper we have presented results that may assist orbital flight dynamics and attitude control 
systems engineers in the design and analysis of INS-equipped spacecraft, by showing how phenomenon 
that occur in terrestrial INS generalizes to spacecraft applications. We showed that the period of the 
Schuler pendulum, the Schuler period in INS error propagation, the LEO orbital period, and the period of 
the oscillatory part of the well-known Hill-Clohessy- Wiltshire equations of LEO satellites' relative motion 
are identical. We have also shown that generalizing the INS error equations results in Hill's equations. 
The satellite relative position in this case corresponds to the INS position error, thus the Schuler frequency 
that is present in the errors of the level channels of terrestrial INS generalizes to an out-of-plane relative 
position oscillation at the orbit frequency. Although these oscillations do not directly enter into the attitude 
errors of a spacecraft-mounted INS, they do affect pointing of the spacecraft relative to the earth, since 
their effect on the navigation states corrupts the spacecraft’s on-board computation of the local level 
reference frame. We also described how INS gyro drift produces an additional Schuler oscillation in the 
attitude error when the INS reference coordinate system is the orbital local-level system. However, when 
the INS reference system is inertial, the attitude exhibits unbounded error growth due to gyro drift. 
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